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laser  beam  propagating  in  a turbulent  medium.  These  results  have  been  used 
to  calculate  the  probability  that,  due  to  beam  •srander.  the  beam  •arill  fail  to  bit 
a receiver  aperture. 
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Short-Term  Spot  Size  and  Beam  Wander 
in  a Turbulent  Medium 


1.  INTRODUCTION 

One  difficulty  in  using  laser  systems  for  communication  through  the  earth’s 
atmosphere  is  that  there  is  a possibility  that  the  beam  may  carrier  from  position 
to  position,  and  may  therefore  completely  miss  the  receiver  in  some  instances. 
Since  the  spot  dances  from  position  to  position  in  fames  of  order  1 millisecond, 
this  would  not  be  a problem  for  systems  using  poises  much  longer  than  1 milli- 
second. bat  voald  be  an  important  consideration  for  systems  using  pulses  short 
compared  vith  a millisecond,  in  this  report  we  shall  therefore  study  the  nature 
of  the  beam  wander  and  derive  results  appropriate  for  application  to  laser  com- 
munications systems. 


2.  DERIVATION  OF  BEAM  - WANDER  FORMULA 

It  is  well  understood  from  purely  physical  considerations  that  eddies  which 
are  large  compared  with  the  diameter  of  a laser  beam  tend  to  deflect  the  beam, 
while  those  smaller  than  the  beam  tend  to  broaden  the  beam  but  do  not  deflect  it 
significantly.  Let  us  consider  the  nature  of  the  received  spot  on  an  aperture  in  a 
turbulent  medium.  If  we  loot  over  very  short  times  we  see  a broadened  spot:  as 
we  look  for  longer  times  we  see  that  the  spot  dances  from  position  to  position. 
Therefore  if  we  average  the  received  intensity  over  very  long  times,  the  total 
broadened  spot  would  consist  or  two  components:  actual  short-term  beam 
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■rhere  eQ  is  the  aperture  field  oi  the  transmitter.  L is  the  path  length  in  turbu- 
lence. kQ  is  the  atavenumber.  and  Mgj  is  the  short-term  rsodulation  transfer 
function  given  by*1  (assuming  D > I where  (Q  is  the  inner  scale  of  the  turbulence) 


M 


s-p  ip.  L)  = exp 


-0.13  *2k  2p5f3 
o 

. 5*3 


(2) 


where  y is  a number  of  order  unity  and  C“  is  the  index  of  refraction  structure 

n 

constant.  In  Eq.  (2)  the  integral  cn  £ from  y (p/D)  to  « indicates  that  the  effect 
on  the  MTF  of  eddies  larger  than  D is  excluded.  For  y ( p/D ) < 1 . Eq.  (2)  can 
be  approximated  by  (assuming  also  that  D < Lq  where  Lq  is  the  outer  scale  size 
of  the  turbulence) 


For  y i p/D)  » 1 . re  can  approximate  by 


(3b) 


(4) 


For  yjjjjof  order  1 . we  have  not  been  able  to  derive  any  approximate  results  for 
Mg— . and  the  fall  expression  of  Eq.  (2)  must  be  used. 

The  long-term  averaged  irradiance  <I>LT  is  given  by  Eq.  (1)  with  re- 
placed by  the  long-term  modulation  transfer  function  M*  is  given  by 

Eq.  (2)  with  y = 0 . and  is 


(5) 


4.  Clifford.  S. . Ochs.  G. . and  Lawrence.  R-  (1974)  Saturation  of  optics  1 scin- 
tillation by  st -oag  tarbuleace.  J.  Opt.  Soc.  An.  S4:148-154. 
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In  general,  the  expressions  for  the  long-  and  short-term  irradiance  most  be 
evaluated  numerically,  as  are  have  done  in  Section  3.  Hen: ever,  there  are  limiting 
cases  where  approximate  analytical  results  can  be  obtained. 

(a)  >5>  I - In  this  cast  Mgr  =*  * *.  and  we  get 


<!>sr  - <I>LT  * -vd>  Jo  f~--' 

0 


“ *VACOUM 

where  we  have  assumed  an  aperture  distribution 


e = exo 
o 


44 


and  have  defined 


a = 


k0°p 


9 4 

koD 
16  L2 


•*6» 


(7) 


Therefore  if  the  coherence  length  pQ  is  much  greater  than  the  a pert  art  diameter, 
the  long-  and  short-term  averaged  irradiances  at  the  receiver  are  the  same  as 
would  be  present  at  the  receiver  if  the  turbulent  medium  were  replaced  by  vacuum. 

fb)  TPojte  « 1 . In  this  case  the  short-  and  long-term  averaged  irradiances 
can  be  approitfra.ied  by 


<Up)\r 


J 


ydyJ^teyi  exp 


S-^I— i- 


y2(l+^2)|.  (8) 


LT 


- 23 


J Kvi  exp 
»» 


-y2(l  + 02lj. 


(9) 


Equations  (8)  and  (91  are  readily  evaluated  numerically, 
to  estimate  the  radius  of  the  beam  by  approximating 


exp 


( i rW 

m [ i/3]| 

I.  my 

1 W 

) [l-0.67(yy)  j 

However,  it  is  possible 


8 


(10) 


2/  2 

bv  exp  ( - v / v ).  rtcrt  v satisfies 
* ■ / ■ o o 


If  xe  use  this  approximation  -re  have 


<1<P>>«rr  - 


57 


exp 


I 

1 4(l  + j32+yJ2)j 


(11) 


F rom  Eg.  (10)  -re  see  that  the  mean  square  beam  radius  is 


O O 

9 

4L“  D~ 

4L“ 

9 9*  4 

k~D~  * 

o 

9 9 9 * 

k^D“v* 

O - o 

If  -re  solve  Eq.  (101  for  ?Q.  -re  find  (for  « 11 


(12) 


(13) 


(14) 


Equation  (14)  gives  the  mean  square  short-term  beam  radius.  The  long-term 
mean  square  beam  radius  is  readily  obtained  from  Eq.(14)  by  sevting  y = 0. 
'Ve  get 


4 L~ 
2 2 
koD“ 


ALT 
o o • 

W 


(15) 


From  Eqs.  (14)  and  (15)  it  is  nor-  possible  to  calculate  the  mean  sruare  beam 
rander  o„. . since  the  mean  square  beam  -rander  is  related  to  and  o^j.  as 
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(16) 


2 2 2 
°lt  “ °sr  + °v: 


If  we  ase  Eqs. (14)  and  (151  in  (IS),  we  get 

Since  we  have  been  assuming  that  (yp  /Dl  « 1 . we  can  approximate  Eq.  (17) 

w 

farther  by 


(17) 


if?  ’ 

,2  3.24  y^  L~ 

W kiDli3  5-3 

o o 


.1/3  .2 


= 4.70  y L (■ 


.1/3  .5*3 


dz'  C2(z')  (z)^ 
n 


(18) 


2 2 2.  2 

Kow  O^j.  = L where  <6  > is  the  _ jean  square  angle  of  wander  of  the  beam 

We  therefore  have  for  >Po^5  « i : 


<♦*» . 


.1/3 


.5/3 


(19) 


We  have  rot  yet  specified  y1^.  but  estimates  obtained  by  comparison  with  data 

appear  to  indicate  ‘hat  0.9  < ylQ  < 1.7 . We  shall  (moose  y1*3  = 1.0  . since  this 

gives  a reasonable  fit  with  the  data  of  Dowling  and  Livingston.5  and  vith  theo- 

2 

ret j cal  results  of  Fried.” 


3.  NUMERICAL  RESULTS 


Equations  (14)  and  (15)  give  the  short-  and  long-term  beam  spread  for 
yo^  « l . However,  we  will  also  often  be  interested  in  the  case  when 
TP^)  — 1 - In  this  case  the  long-  and  short-term  ir radiances  must  be  evaluated 
numerically.  If  we  assume  T.  c tperture  field  is  given  by  Eq.  (7).  then  we  must 
evaluate 


<I(p)>ST  = 202  J ydyJo(ay)  expj-y* 


u+p  i 


(20) 


5.  Dowling.  J.  and  Livingston.  P.  (1973)  Behavior  of  focussed  beams  in  atmo- 
spheric turbulence:  Measurements  and  comments  on  the  theory.  J.OpL 
Soc.  Am.  63:846-858.  
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Equations  (20)  and  (21)  have  been  evaluated  numerically  to  determine  the  long- 
term and  short-term  beamwidth  as  a function  of  yp0jD.  for  various  values  of 
when  r=l.  These  results  are  shown  in  Figure  2.  From  this  figure  we  can 
readtlv  deter. nine  the  mean  square  beam  wander  since 


random  variable,^  we  can  nsv  calculate 
point  aperture.  We  have  for  the  prob- 


12 */ 

The  probability  that  the  beam  will  miss  (by  miss  we  mean  that  the  amplitude  at 
the  receiver  will  be  reduced  by  more  than  e *)  a point  aperture  is  equal  to  the 
probability  that  the  beam  wander  exceeds  the  short-term  beam  spread.  That  is 


Since  the  beam  wander  is  a gaussian 
the  probability  that  the  beam  will  miss  a 
ability  of  the  beam  wander 


P <PW>  = 


2to, 


exp 


W 


2 

Pw 

*»w 


ao 


Upon  using  Eq.{23)  in  (25)  we  have 


This  result  io  plotted  in  Figure  3.  For  th«?  case  of  a receiving  aperture  of 
radius  AR.  we  replace  Ogy  h>  in  So.  (SI.  In  the  next  section  we 

will  da  an  example  to  illustrate  the  use  of  Figures  2 and  3 for  atmospheric 
propagation. 


4.  ATfLlCVnON  TO  THE  EAKllfS  ATMOSfUERE 

An  important  application  of  our  results  is  to  a laser  communications  link 
between  the  earth  and  an  orbiting  satellite.  Let  us  suppose  the  satellite  is  at  an 
altitude  hs  and  at  an  angle  6 relative  to  the  normal  to  the  earth  at  the  transmitter, 
as  shown  in  F;gure  4.  For  this  case  we  may  write 


12 


2 


(Sst/sLt) 

Figure  3.  Probability  That  a Laser  Beam  Will  Miss  a 
Point  Receiving  Aperture,  as  a Function  of  the  Ratio  of 
the  Short-Term  to  Long-Term  Spot  Size 


(h  -*•  a)  dh 

Vh2  2ah  a"  cos~0 


Tfhere  a is  the  radius  of  the  earth,  if  ve  assume  h « a and 


(28) 


d < cos 


-1 


\ll 


V a / 


,1/21 


we  can  approximate  Eq.  (27)  bjr 


dr  = sec  6 dh  . 


(29) 


and  we  nary  therefore  evaluate  pQ.  for  the  earth-to- satellite  path  as 


-5*3 


1.46  k2  sec  6 r-hs 


£ * dh  C2(h)  (M5* 


(30) 


2 

A common!*.-  used  distribution  for  C is 

n 


(31) 


Cn  " Co  b'>’  [-h'»o] 


so  that  for  h » h 
s c 


-5/3  _ 2 -p/3  v 2.S&-V.-5/3 

p = 1.46  k sec  c - S’ i C fa  b 
o o »3  * o o s 


(32) 


where  T is  the  gamma  function.  If  we  take  - 10  *3.  v - 1/3  and  hQ  = 1000 
meters,  we  get 


3/5 


pQ  = hs  J8.C3  X 10  9 k “ sec  fij  meters. 


(33) 


For  8 = 30*  and  kQ  = 0.95  X 107  m 1 (0.633  pin  light),  >e  have 


p„  = 2.63  X 10  h meters, 
o s 


(34) 


Now  suppose  the  transmitter  diameter  D is  1 m a*ri  the  satellite  is  at  an 
altitude  hs  of  200  km.  Then  p ^ = 52.5  and  p = 10  . Consequently,  from 

Figure  2 we  have-  <r«- r^LT=  * aat^  ' 1-0  33  ngore  3 we  see  that  the  probability 
that  the  beam  (over  times  of  order  D/v)  will  miss  a point  receiver  on  the  satel- 
lite is  =*  0.  Of  coarse  a s 8 •*  50*.  P-t^cc  *iil  co  longer  be  nearly  equal 


to  zero. 
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